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Abstract. We have studied bounds on the neutrino mass using new data from the 
WMAP 3 year data, the Sloan Digital Sky Survey measurement of the baryon acoustic 
peak, the Type la supernovae from SNLS, and the Lyman-a forest. We find that 
even in the most general models with a running spectral index where the number of 
neutrinos and the dark energy equation of state are allowed to vary, the 95% C.L. 
bound on the sum of neutrino masses is m v — 0-62 eV (95% C.L.), a bound which 
we believe to be robust. In the more often used constrained analysis with N v = 3, 
w = — 1, and a s — 0, we find a bound of 0.48 eV without using the Lyman-a data. 
If Lyman-a data is used, the bound shrinks to J2 m f — 0-2 — 0.4 eV (95% C.L.), 
depending strongly on the Lyman-a analysis used. 
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1. Introduction 

In the past few years a new standard model of cosmology has been established in which 
most of the energy density of the Universe is made up of a component with negative 
pressure, generically referred to as dark energy. The simplest form of dark energy is the 
cosmological constant, A, which obeys Pa = — Pa- This model provides an amazingly 
good fit to all observational data with relatively few free parameters and has allowed 
for stringent constraints on the basic cosmological parameters. 

The precision of the data is now at a level where observations of the cosmic 
microwave background (CMB), the large scale structure (LSS) of galaxies, and Type 
la supernovae (SNIa) can be used to probe important aspects of particle physics such as 
neutrino properties. Conversely, cosmology is now also at a level where unknowns from 
the particle physics side can significantly bias estimates of cosmological parameters. 

The combination of all currently available data from neutrino oscillation 
experiments suggests two important mass differences in the neutrino mass hierarchy. 
The solar mass difference of Am 2 2 ~ 8 x 1CT 5 eV 2 and the atmospheric mass 
difference Am| 3 ~ 2.2 x 1CT 3 eV 2 [1]. In the simplest case where neutrino masses are 
hierarchical these results suggest that mi ~ 0, rri2 ~ Am so i ar , and m 3 ~ Am atm ospheric • 
If the hierarchy is inverted one instead finds m 3 ~ 0, m 2 ~ Am atmosp h e ric , and 
mi ~ Am atm0 sphcric- However, it is also possible that neutrino masses are degenerate, 
mi ~ m 2 ~ m 3 ^> Am atmosphcric . Since oscillation probabilities depend only on squared 
mass differences, Am 2 , such experiments have no sensitivity to the absolute value of 
neutrino masses, and if the masses are degenerate, oscillation experiments are not useful 
for determining the absolute mass scale. 

Instead, it is better to rely on kinematical probes of the neutrino mass. Using 
observations of the CMB and the LSS of galaxies it has been possible to constrain 
masses of standard model neutrinos. The bound can be derived because massive 
neutrinos contribute to the cosmological matter density, but they become non-relativistic 
so late that any perturbation in neutrinos up to scales around the causal horizon at 
matter-radiation equality is erased, i.e. the kinematics of the neutrino mass influences 
the growth of structure in the Universe. Quantitatively, neutrino free streaming 
leads to a suppression of fluctuations on small scales relative to large by roughly 
AP/P ~ — &VL u /VL m [2]. The density in neutrinos is related to the number of massive 
neutrinos and the neutrino mass by 

n , 2 _ TjUv _ N v m u 

v ~93.2eV~93.2eV K) 
where h is the Hubble parameter in units of 100 km s~ x Mpc -1 and all neutrinos 
are assumed to have the same mass. Such an effect would be clearly visible in LSS 
measurements, provided that the neutrino mass is sufficiently large, and a likelihood 
analysis based on the standard ACDM model with neutrino mass as an added parameter 
in general provides a bound for the sum of neutrino masses of roughly X) m ^^0.5 — 1 
eV, depending on exactly which data is used [4-14] . 
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This should be compared to the present laboratory bound from 3 H beta decay 
found in the Mainz experiment, m Ue = ^^mf) 1 ^ 2 < 2.3 eV [15]. It should also be 
contrasted to the claimed signal for neutrinoless double beta decay in the Heidelberg- 
Moscow experiment [16-18], which would indicate a value of 0.1-0.9 eV for the relevant 



combination of mass eigenstates, m ee = J2j U^m^ 



Some papers claim that the 



cosmo logical neutrino mass bound is already incompatible with this measurement. 

However, those claims are based on a relatively limited parameter space. Using 
a much more complicated model with a non-power law primordial power spectrum it 
is possible to accomodate large neutrino masses, provided that SNIa data is discarded 
[40] . It was recently shown that there is a strong degeneracy between neutrino mass and 
the equation of state of the dark energy component when CMB, LSS, and SNIa data is 
considered. The reason is that when the neutrino mass is increased, the matter density 
must be increased accordingly in order not to conflict with LSS data. This is normally 
excluded when the dark energy is in the form of a cosmological constant. However, 
when the equation of state parameter of the dark energy fluid, w, is taken to be a free 
parameter, an increase in the matter density can be compensated by a decrease in w to 
more negative values. 

Here we study how the degeneracy can be broken by the addition of information 
provided by the measurement of the baryon acoustic beak in the Sloan Digital Sky 
Survey data and the distances to Type la supernovae provided by the SNLS data. We 
also study constraints derived from combining all available cosmological data, including 
the measurement of the small scale matter power spectrum from the Lyman-a forest. 

In the next section we discuss the cosmological data used, in section 3 numerical 
results from the likelihood analysis are provided, and finally section 4 contains a 
discussion. 



2. Cosmological data 

In order to probe the neutrino mass we 
CMB, LSS and SNIa data. 



have used some of the most recent precision 



2.1. Type la supernovae. 

Observations of SNIa over a wide redshift range to determine cosmological distances 
is perhaps the most direct way to probe the energy content of the universe [19] 
and have lead to a major paradigm shift in cosmology [20-25]. While extremely 
succesful at showing the need for dark energy to explain the derived distances, these 
data sets are plagued by systematic uncertainties (e.g. dust extinction corrections, 
K- correct ions, calibration uncertainties, non-la contamination, Malmquist bias, weak 
lensing uncertainties, etc.) rendering them non-optimal for precision tests on w. 
However, the situation has improved since the start of a dedicated experiment, the 
SuperNova Legacy Survey (SNLS) at CFHT. Thanks to the multi-band, rolling search 
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technique, extensive spectroscopic follow-up at the largest ground based telescopes and 
careful calibration, this data-set is arguably the best high- ,2 SNIa compilation to date, 
indicated by the very tight scatter around the best fit in the Hubble diagram and the 
careful estimate of systematic uncertainties by Astier et al., (2006) [26]. The first year 
data of the SNLS collaboration includes 71 high redshift SNIa in the redshift range 
z = [0.2, 1] and 44 low redshift SNIa compiled from the literature but analysed in the 
same manner as the high-z sample. We thus make only use of this new SNIa data set 
to minimize the effects from systematic uncertainties in our analysis. 

It should be noted that combining the results in Ref. [26] with the Gold sample of 
Riess et al. (2004) [25] as done in, e.g. Refs. [3] and [62] is not straight forward. First, 
the two data-sets use the same low-redshift samples to anchor the Hubble diagram and 
are thus not independent. Second, there are systematic differences in the way the SNe 
are analysed, e.g. with respect to reddening corrections. In a combined analysis, this 
would have to be adressed. 

2.2. Large Scale Structure (LSS). 

Any large scale structure survey measures the correlation function between galaxies. 
In the linear regime where fluctuations are Gaussian the fluctuations can be described 
by the galaxy-galaxy power spectrum alone, P(k) = \5k, gg \ 2 - In general, the galaxy- 
galaxy power spectrum is related to the matter power spectrum via a bias parameter, 
b 2 = P gg /P m . In the linear regime, the bias parameter is approximately constant, so up 
to a normalization constant, P gg does measure the matter power spectrum. 

At present there are two large galaxy surveys of comparable size, the Sloan Digital 
Sky Survey (SDSS) [27, 28] and the 2dFGRS (2 degree Field Galaxy Redshift Survey) 
[29]. Once the SDSS is completed in 2006 it will be significantly larger and more accurate 
than the 2dFGRS. In the present analysis we use data from both surveys. In the data 
analysis we use only data points on scales larger than k = 0.15/i/Mpc in order to avoid 
problems with non-linearity. 

2.3. Baryon acoustic oscillations (BAO). 

The acoustic oscillations at the time of CMB decoupling should be imprinted also on 
the low-redshift clustering of matter and manifest themselves as a single peak in the 
galaxy correlation function at ~ 100/i" 1 Mpc separation. Because of the large scale and 
small amplitude of the peak, surveys of very large volumes are necessary in order to 
detect the effect. 

A power spectrum analysis of the final 2dFGRS data shows deviations from a 
smooth curve at the scales expected for the acoustic oscillations. The signature is much 
smaller than the corresponding acoustic oscillations in the CMB but can be used to 
reject the case of no baryons at 99% C.L. [30]. 

The SDSS luminous red galaxy (LRG) sample contains 46,748 galaxies with 
spectroscopic redshifts 0.16 < z < 0.47 over 3816 square degrees. Even though the 
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number of galaxies is less than the 2dFGRS or the main SDSS samples, the large survey 
volume (0.72/i~ 3 Gpc 3 ) makes the LRG sample better suited for the study of structure 
on the largest scales. The LRG correlation function shows a significant bump at the 
expected scale of ~ 150 Mpc which, combined with the detection of acoustic oscillations 
in the 2dFGRS power spectrum, confirms our picture of LSS formation between the 
epoch of CMB decoupling and the present. 

For a given cosmology, we can predict the correlation function (up to an amplitude 
factor which is marginalized over) and compare with the observed LRG data. The 
observed position of the peak will depend on the physical scale of the clustering and 
the distance relation used in converting the observed angular positions and redshifts to 
positions in physical space. The characteristic physical scale of the acoustic oscillations 
is given by the sound horizon at the time of CMB decoupling and depends most strongly 
on the combination Vt m h 2 . The conversion between positions in angular and redshift 
space to positions in physical space will cause the observed correlation scale to depend 
on the distance combination 

nl/3 



D v (z) = 



D M {zf ' 



(2) 



H{z) 

where H(z) is the Hubble parameter and Dm(z) is the comoving angular diameter 
distance. 

Our approach to fit the data has been to calculate the matter power spectrum for 
a given model, then Fourier transforming it to obtain the 2-point correlation function, 
£(r). This correlation function has been fitted to the SDSS data using the full covariance 
matrix given by [31]. 

In terms of the simple parametrization provided by [31] in terms of the parameter 



A = D v (z)^l, (3) 

zc 

we find that the SDSS constraint can approximately be written as 

/ fl \ -0-35 

A = 0.469 (—J (1 + 0.94/,) ±0.017, (4) 
where f v = Vt u /Vt m . 



2.4- The Lyman-a forest. 

Measurements of the flux power spectrum of the Lyman-a forest has been used to 
measure the matter power spectrum on small scales at large redshift. By far the 
largest sample of spectra comes from the SDSS survey. In Ref. [33] this data was 
carefully analyzed and used to constrain the linear matter power spectrum. The derived 
amplitude is A 2 (k = 0.009 kms -1 , z = 3) = 0.452io!o6 an d the effective spectral index 
is n e ff = — 2.3211q;o5. The result has been derived using a very elaborate model for 
the local intergalactic medium, including full N-body simulations. It has been shown 
that using the Lyman-a data does strengthen the bound on neutrino mass significantly. 
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However the question remains as to the level of systematic uncertainty in the result. 
The same data has been reanalyzed by Seljak et al. [62] and Viel et al. [63-65], with 
somewhat different results. The normalization found in Refs. [63-65] is lower than that 
found by Ref. [33] which has significant influence on the inferred neutrino mass bound. 

Our Lyman-a analysis uses only the effective parameters A 2 (A;) and n s . The 
resulting bound on m u using these parameters was shown in Ref. [12] to be very similar 
to those obtained in a full analysis. Furthermore we stress that using these parameters 
makes it very simple to test the influence on different assumptions about the Lyman-a 
data on the neutrino mass bound. 

2.5. Cosmic Microwave Background (CMB). 

The temperature fluctuations are conveniently described in terms of the spherical 
harmonics power spectrum Ct,i = (K m | 2 ), where -^(#,0) = Y,i m o-imXimiO, </>)• Since 
Thomson scattering polarizes light, there are also power spectra coming from the 
polarization. The polarization can be divided into a curl-free (E) and a curl (B) 
component, yielding four independent power spectra: Ct,i, Ce,i, Cb,i, and the T-E 
cross-correlation Cte,i- 

The WMAP experiment has reported data on Ct,i and , Cee,u an d Cte,i as 
described in Refs. [3, 34, 35] . We have performed our likelihood analysis using the 
prescription given by the WMAP collaboration [3,34,35]. 

We furthermore use the newly published results from the Boomerang experiment 
[37-39] which has measured significantly smaller scales than WMAP. 

3. Likelihood analysis 

Using the presently available precision data we have performed a likelihood analysis for 
the neutrino mass. 

As our framework we choose a flat dark energy model with the following free 
parameters: Q m , the matter density, the curvature parameter, f^, the baryon density, 
w, the dark energy equation of state, H , the Hubble parameter, n s , the spectral index 
of the initial power spectrum, a s , the running of the primordial spectral index, N v , the 
effective number of neutrino species, and r, the optical depth to reionization. Finally, the 
normalization, Q, of the CMB data, and the bias parameter b are used as free parameters. 
The dark energy density is given by the flatness condition Qde = 1 — fi m — fij,. Including 
the neutrino mass our benchmark model has 11 free parameters. We also test more 
restricted parameter spaces in order to probe parameter degeneracies. 

The priors we use are given in Table I. The prior on the Hubble constant 
comes from the HST Hubble key project value of h = 0.72 ± 0.08 [36], where 
ho = i/o/lOOkms^Mpc" 1 . 

When calculating constraints, the likelihood function is found by minimizing \ 2 
over all parameters not appearing in the fit (i.e. over all other parameters than m„). 
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Table 1. The different priors on parameters used in the likelihood analysis. 

3.1. Results 

In Fig. 1 we show the one dimensional likelihood function for the neutrino mass for 
various different data sets, using the full 11-dimensional parameter space. 

By far the most conservative is for the case which includes only CMB, LSS, and 
SN-Ia data. Because the parameter space is larger than the one used in Ref. [41], the 
constraint of 1.72 eV is, however, stronger because the addition of new WMAP-3 data 
breaks some of the degeneracy f. It should be noted though that the bound is weaker 
than the one shown in Fig. 18 of Ref. [3]. The main reason for this is most likely that 
N v and a s is allowed to vary. 

When the BAO data is added it has the effect of essentially breaking the degeneracy 
between m v and w. The reason is that the BAO measurement is almost orthogonal to 
the SNIa measurement in the [Q m , iu]-plane. For this case the bound shrinks to 0.62 eV, 
a factor of three improvement. 

We also test the inclusion of SDSS Lyman-a data in the fit. When BAO data is not 
included the Lyman-a data itself does break some of the degeneracy, to a level where 
the formal bound is 0.83 eV. Including both BAO and Lyman-a data does lead to a 
significant improvement, the combined bound being 0.49 eV at 2a. 

It is interesting to see how this bound compares with the value obtained in the 
standard ACDM+m^ parameter space (which has 8 parameters). Including the BAO 
data we find an upper bound of 0.48 eV, and with the SDSS Lyman-a data we find 
0.35 eV. Combining all data leads to a bound of 0.27 eV. It should be noted that this 
bound is somewhat weaker than what is found in Ref. [62] , presumably because of the 

| Note, however, that if only WMAP-3 data is used, the inferred bound of 2 eV [42] is almost identical 
to what is found from WMAP-1. This means that if used alone WMAP-3 does not add significant new 
information 
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Data m v (95% C.L.) 
1: CMB, LSS, SNIa 1.72 eV 

2: CMB, LSS, SNIa, BAO 0.62 eV 

3: CMB, LSS, SNIa, Ly-a 0.83 eV 

4: CMB, LSS, SNIa, BAO, Ly-a 0.49 eV 



Table 2. Best fit x 2 values for the four different analyses presented in Fig. 1, in all 
cases based on the full 11-dimcnsional parameter space. 



Data m v (95% C.L.) 
1: CMB, LSS, SNIa 0.70 eV 

2: CMB, LSS, SNIa, BAO 0.48 eV 

3: CMB, LSS, SNIa, Ly-a 0.35 eV 

4: CMB, LSS, SNIa, BAO, Ly-a 0.27 eV 



Table 3. Best fit % 2 values for the three different analyses presented in Fig. 2, in all 
cases based on the restricted 8-dimensional parameter space with N v = 3, w = — 1, 
and a s = 0. 



new treatment of Lyman-a in Ref. [62]. 

The strength of the bound including Lyman-a therefore depends crucially on the 
assumed uncertainty in the measurement of A 2 (k). We have not used the most recent 
analysis of the SDSS data by Seljak et al. [62], but in order to compare roughly with 
their result we have added Lyman-a data that mimics what is shown in their Fig. 1. 
For this model we find a bound of 0.20 eV, in good agreement with their result of 0.17 
eV. 

However, we also note that using the Lyman-a result from Viel et al. [63-65] which, 
when combined with WMAP-3, has a best-fit normalization (note that this is mainly 
due to the fact that the high resolution Lyman-a data used in Refs. [63-65] has larger 
error bars) about 2a lower than in Ref. [62], would lead to a different bound on J2 m v 
To test this we have also run a likelihood analysis using the Viel et al. data on A 2 (k) and 
found a bound of 0.40 eV. This is the limit we inferred by using the derived cosmological 
parameters from the analysis of the SDSS data of [63-65], which consists of a Taylor 
expansion of the flux power around a best fit model (SDSS-d in their tables). With this 
analysis the bound from Lyman-a adds relatively little information to that obtained 
from BAO data. We show the \ 2 curves for this analysis in Fig. 3. 

This result leads to the inevitable conclusion that while the Lyman-a data holds 
the potential to provide crucial information on the neutrino mass, there seem to be 
unresolved systematic issues related to the conversion of the measured flux power 
spectrum to A 2 (A;). Any bound derived from the use of the Lyman-a power spectrum 
should therefore be treated with some caution. This is especially true since in the Seljak 
et al. [62] analysis there is a tension between the WMAP and Lyman-a normalization on 
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Figure 1. The value of A% 2 as a function of ^rn^, for various different data sets 
used for the full 1 1-dimensional parameter space. The curves are identical to the cases 
described in Table 2: The full curve is case 1, the long-dashed is case 2, the dotted is 
case 3, and the dashed is case 4. 

small scales of more than 2a. While this could be statistical in nature, the difference with 
respect to the Viel et al. [65] analysis of the same data suggests a possible systematic 
origin. 

4. Discussion 

We have calculated the bound on the sum of light neutrino masses from a combination 
of the most recent cosmological data. If only CMB, LSS, and SNIa data is used, we find 
a strong degeneracy between Y^ m u, N v , and w which severely limits the ability of this 
data to constrain the neutrino mass. However, once data from the SDSS measurement 
of the baryon acoustic peak is included, this degeneracy is broken because it measures 
Q m and w very precisely. The derived bound and best fit model is compatible with 
what can be derived from observations of the Lyman-a forest, but is likely much 
less affected by systematics. The bound is J2 m u ^0.6 eV even for a very general 11 
parameter cosmological model. If data from the Lyman-a forest is added, the bound is 
J2 m u < 0.2 — 0.4 eV (95% C.L.), depending on the specific Lyman-a analysis used. 

Beyond SDSS, future galaxy redshift surveys will achieve a larger effective volume 
and go to higher redshifts (see for instance Ref. [43]). The prospects for them to constrain 
m v have been studied in e.g. Refs. [44,45]. With the help of the BAO measurement, the 
angular diameter distance and the Hubble parameter H(z) can be measured to percent 
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Figure 2. The value of A% 2 as a function of m v for various different data sets used 
for the restricted 8-dimensional parameter space with N v = 3, w = —1, and a s = 0. 
The curves are identical to the cases described in Table 3: The full curve is case 1, the 
long-dashed is case 2, the dotted is case 3, and the dashed is case 4. 




Figure 3. The value of A% 2 as a function of ^ m„ for the restricted 8-dimensional 
parameter space with N v = 3, w = — 1, and a s — for different assumptions about the 
Lyman-a data. The full line is for the Ref. [33] data, the dashed is for the approximate 
analysis with the data in Ref. [62] , and the long-dashed is for the data in Ref. [65] . 
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level. This enables w to be determined within ~ 10% [46-48], breaking the m v — w 
degeneracy. 

Another powerful probe in the future is weak gravitational lensing. It traces 
directly the mass distribution in a wide range of scales, and thus does not suffer from 
the light-to-mass bias in large scale structure surveys, while still being sensitive to 
effects from neutrino early free-streaming. Cosmic shear measurement with tomographic 
redshift binning of sources galaxies is particularly effective in constraining dark energy 
parameters. The degeneracy between m v and w can then be broken in a similar way as 
BAO does. Systematics herein arise from photometric redshift uncertainties and shear 
calibration errors, which are expected to be under control to the required accuracy. 
The major uncertainty comes from estimating the nonlinear part of the matter power 
spectrum [49,50]. It is found that future ground- and space-based surveys such as 
CFTHLS [55], SNAP [56] and LSST [57], combined with future CMB measurement 
from the Planck Surveyor [58], can constrain neutrino mass to a(J2 m u) = 0.025 — 0.1 
eV [51,52]. 

CMB photons from the last scattering surface at z 1100 are also deflected by 
the large scale structure at z <>3. Extracting the weak lensing information encoded in 
the CMB signal will significantly enhance the sensitivity of CMB experiments to small 
neutrino mass [53,54]. The forecast error o~(^2m u ) obtained in Ref. [53] is ~ 0.15 eV 
for Planck and SAMPAN [59], and 0.035 eV for the future Inflation Probe project [60]. 
To summarize, we expect near future cosmological observations to pin down the total 
neutrino mass to a precision better than 0.1 eV, the level of the inverted hierarchy. 
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